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ABSTRACT 

Context. Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples 
include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where 
all of the photons energetic enough to ionize hydrogen have been absorbed. 

Aims. To determine the physical properties of the PDRs associated with the star forming regions IRAS 23133-1-6050 and S 106 and 
present them in the context of other Galactic PDRs associated with massive star forming regions. 

Methods. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 pm spectrum of each object 
from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements 
(and combinations thereof) are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared 
to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. 

Results. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very 
dense, highly UV irradiated phase (n ~ 10® cm^'^, Go ~ 10®) interspersed within a lower density, weaker radiation field phase (n ~ 10“* 
cm“®, Go ~ 10“*). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak 
of the '^CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical 
modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the *^CO 
ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component 
for the hightest J lines (J ^ 20) is required for S 106. 

Key words. ISM: general; ISM: Molecules; photon-dominated region (PDR); Infrared: ISM; Stars: massive; Stars: formation 


1. Introduction 


Understanding photodissociation regions (PDRs; also known as 
photon dominated regions; e.g. iTielens & Hollenbachl 1 198^ 
ISternberg & Dalgarnoll98^ is an essential part of comprehend¬ 
ing the physics and chemistry of the interstellar medium (ISM). 
The classical definition of a PDR is the region in which the en¬ 
ergy balance is dominated by far-UV photons with energies less 
than that required to ionize hydrogen (13.6 eV). In practice this 
definition means that classical PDRs form the ‘skin’ of the ion¬ 
ized H II regions and represent the buffer zone between the fully 
ionized H ii region material and molecular clouds, typically this 
means that they occur at Ay < 5 with respect to the cloud sur¬ 
face. In a broader sense, a large fraction of all neutral ISM mate¬ 
rial exists under these conditions, underlining the importance of 
understanding the physics and chemistry of these regions. 

The prototypical PDR, th e Orion (Bar) region, has 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 
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has been largely due to two factors: a) its close proximity (450 
PC. iGenzel & StutzOll989l) and b) the fact that the Orion Bar 
is viewed almost exactly edge on. This combination allows for 
the relatively easy separation of the emission from each region 
across the transition from ionized to PDR materi al and subse¬ 
quent comparison with detailed PDR models (e.g. . ITielens et al.l 
Il993l) . These models are highly sophisticated codes that include 
a detailed chemical network as well as descriptions of relevant 
heating and cooling processes. As a function of the incident 
radiation field strength (usually denoted GcQ) and density (or 
thermal pressure) the PDR models calculate the thermal equi¬ 
librium temperature and dominant atomic and molecular species 
in chemical balance. The models also predict the dominant emis¬ 
sion line strengths that can be compared with observations. From 
this data the PDR models synthetic spectra are constructed and 
compared to those observed. It is now generally agreed that most 
PDRs are clumpy, that is, they take the form of dense clumps 
(n ^ 10^ cm ®) e mbedded within an interclump medium 10-100 
times less dense (iBurton et al.lll99()bt iMeixner & Tielensll 19931: 
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iRollig et'^l2006h . Such two phase models can reproduce many 
of the observed atomic and molecular emission from PDRs. 

Measurements of PDR diagnostics have been sparse in pre¬ 
vious decades because features diagnostic of PDR physical con¬ 
ditions predominantly emit in the far-infrared (FIR) spectral 
regime at wavelengths between ~30 and 1000 pm. This wave¬ 
length region is notoriously difficult to observe, as much of 
this range requires above atmosphere observations, using space- 
based observatories such as ISO, Spitzer and Herschel, bal¬ 
loons such as BICE (iNakagawa et alJ[l998h or aircraft like the 
KAO (Kuiper Airborne Observatory) and SOFIA (Stratospheric 
Observatory for Infrared Astronomy) to cover the shortest 
wavelengths and ground based sub-mm observatories to cover 
the longer wavelength emissions. In recent years this has be- 
come much easier though, as the Herschel space observatory 
(iPilbratt et al.ll2010h was able to collect full spectra in the PDR- 
significant range. 

In this paper we expand the sample of PDRs ob¬ 
served by Herschel, using the PACS and SP IRE instru¬ 
ments (iPoglitsch et al.l 1201 Ot iGriffin et al.l 1201 Oh . to include 
two Galactic PDRs generated by massive star formation 
(IRAS 23133-1-6050 and S 106) in order to understand the phys¬ 
ical conditions and compare their properties to those of other 
Galactic PDRs, including the prototypical Orion Bar. We ob¬ 
tained spectra of the important atomic PDR cooling lines ([O i] 
63, 145 pm, [C ii] 158 pm, and [C i] 609 pm) as well as the 
CO and '^CO ladders from mid- to high-J. These observations 
allow us to probe various depths into the respective PDRs, with 
the atomic lines providing constrains on the warm, dense sur¬ 
face layers at low Ay (e.g. T = 100-1000 K, n = 10"^-10^ cm“^. 
Ay < 4) and the CO ladders probing both the warm, low Ay 
regions, as well as the deeper, cooler regions at higher Ay (e.g. 
T = 10-100 K, n = 10"^-10^ cm“^. Ay > 4 ). The nature of our 
targets and prior knowledge of them is outlined in Section |2l A 
description of the observations, data reduction procedures and 
the measurement of the spectra is given in Section [3] The re¬ 
sults of this work, starting with the PDR diagnostics including 
the molecular tracers such as '^CO are discussed in Section |4] 
in terms of the physical conditions of these PDRs. Subsequently 
we compare these observations to the outputs of a sophisticated 
PDR model in Section|5] Discussion of these results is presented 
in terms of the geometry of the PDRs and the derived condi¬ 
tions with those of other Galactic PDRs in Section|6] Finally, we 
present our conclusions and a short summary in Section |7] 


2. IRAS 23133+6050 and S 106 


Both IRAS 23133+6050 and S 106 are well-studied regions of 
star formation. In this section we review the pertinent litera¬ 
ture in order to guide our later discussion. The properties drawn 
from the literature are summarized in Table[I] GLIMP SE survey 
images (iBeniamin et al.ll^03HChurchwelljr^ 2009b, obtained 


using the Spitzer/IRAC instrument (iFazio et al. 
source are shown in Fig. [1] 


2004b . of each 


* Go is defined in units of the Habing field after iHabind lll968h . in 
which the local interstellar UV radiation field integrated between 6 and 

13.6 eV was found to have an average strength of 1.2 x 10"** ergs cm“^ 
s * sr“*. The mean interstellar radiation fields dPrainel 1978b . is greater 
than the Habing field by a factor of 1.7. 


2.1. IRAS 23133+6050 

IRAS 23133+605CB appears to be a classical example of 
an is olated star formin g regio n within a molecular cloud 
(e.g., iKurtz et ^ 11994 1 19991 show it to be an isolated 
ultra-compact H ii region). From an astrometric perspec¬ 
tive, there has been some debate re garding the distan ce of 
IRAS 23133+6050. Early estimates bv iHumphrevsl (1 1 978h gave 
the distance to the Cassiopeia OB2 star for ming region, o f which 
IRAS 23133+6050 is part, at 2.63 kpc. iPeeters et alJ (l2002b 
found a distance of 5.5 kpc using kinematical methods and ISO- 
SWS spectro s copy, which agrees with an earlier calculation by 
iFich & Blitzl (1 1984b who found d = 5 .2 + 0.9 kpc. Howeyer 
other authors, e.g.. iLebron et al.l (1200lb . adopt a smaller yalue 
of 3.1 kpc, closer to initial estimates, based on s pectrophotmet- 
ric obs eryations of associated stars as deriyed by lBrand & Blitzl 
(Il993b. We ado p t the h igher yalue throughout this paper. 

iKurtz et ^ d 1999b showed that the radio emission from 
IRAS 23133+6050 is consistent with the H ii region being pow¬ 
ered by a late O star of spectral type 09.5. This star powers a 
cometary ultra compact H ii re gion with a diameter of 0.2 pc 
(iMartln-Hernandez et alJ f2003b based on their calculated dis¬ 
tance o£_5;5_^2cMJsing_JTOob^yations of the MIR [S III] 
lines IMartln-Hernandez et al.l (l2003b found the electron density 
of the H II region to be about 2 x 10^ cm“^. Comparatiyely 
few studies haye discussed the PDR / molecular cloud enye- 
lope of the IRAS 23133+6050 region. Assuming pressure equi¬ 
librium, the expected density contrast between a typical H ii 
region and its PDR is expected to be around 30; which leads 
to an expected PDR density of the order of 6 x 10"^ cm“^. 
IRAS 23133+6050 was included in a sample of UC-H ii regions 
obserye d in the sub-mm by SCU BA, proyiding fluxes at 450 and 
850 fim (iThompson et al.ll2006b . 


2.2. S 106 


S 106 (Sh 2-106; Ish^ iii Il959h is an extremely complex 
source, consisting of multiple components and IR sources. 
S 106 has been foun d to be at a distance of 1. 2-1.8 kpc 
( Schneider et al.ll2007b . we adopt 1.5 kpc following iNoel et ^ 
(l200.5b . S 106 is be heyed to be power ed mainly by a single O- 
type star (S106-IR; lRicher et al.l ll99 ^. which is surround ed by 
a yery small optically thick disk ( Hoare & Muxlowlll996b . The 


double lobed structure of the H ii region is created by the shad¬ 
ows of this disk. 

The PDR and molecular cloud components of S 106 haye 
been the focus of prior studies (see iHodapp & Schneideill2008l 
for a thorough reyiew). The oyerall structure of the S 106 sys¬ 
tem is yery similar to bipolar nebulae formed by a combina¬ 
tion of stellar outflows and a circumstellar disk. That the central 
star i s undergoing significant mass loss ha s long been known 
(e.g., iDrew et ^119931: iGibb & Hoarell2()(y^ iParkin et al.ll2009h 
which is consistent with the idea that the large ionized lobes were 
sculpted in this way. Howeyer the origin of the central dark lane 
has been more controyersial. 

Early models of the S 106 system inyoked a large optically 
thick torus of dense material as the explanation for the double 
lobed appearance of the system and central dark lane. Higher 


^ We note that while the designation IRAS 23133+6050 technically 
refers to the exciting star, throughout this paper we will also refer to its 
H II region / PDR complex by this name as appears to be the common 
practice in the literature. The H ii region c an be found in SIM BAD 
as ‘[KCW94] 111.612+0.374’ reflecting the IKurtz et all dl^ radio 
detection of the region. 
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Table 1. Source Characteristics 



S 106 

Reference 

IRAS 23133+6050 

Reference 

Astrometric: 





Right Ascension [h:m:s] 

20:27:26.8 


23:15:33.1 


Declination [°: "] 

+37:22:49 


+61:07:18 


Angular Size of H ii region ["] 

~ 120 X 60 

1 

9 

4 

Distance [pc] 

1500 

2,3 

5500 

4 

Physical Radius of H ii region [pc] 

0.44“ 

2 

0.1 

4 

Physical Radius of Herschel Beam [pc] 

0.15 


0.49 


Central Stars: 





Name 

S106-IR 

5 

IRAS 23133+6050 


Spectral Type 

08 

6 

08.5 

9 


BO-08 

7 

09.5 

10 


07-09 

8 

08.5 

11 

ISM Properties: 





Electron density [cm“^] 

1.5 -6 X 10^ 

12 

2 X 10® 

2 




8.6 X 10® 

13 

PDR density [cm“^] 

1 -3.5 X 10^ 

1 




0.3-5 X 10® 

14 




~ 10® 

4 




Here we have adopted 60" as an average. 


Re ferences. tlllSi mon et a f](l2012h : f2ilNoel et alj(l2005h: GllSchneider et al! 1200'^ : t41lMartm-Hemandez et alj(l2003 1: (5) Richer et '5](ll993h: 
l lvan den Ancker et all l l20Q(Hh fTi lEiroa et alldlQT^ Si lGehrz et al.l(Il982h: (9ilKurtz et alHl999ll : (lOi lKurtz et al.l(11994) : G li lHanson et al.l 


11200 


2002h : fl2i llsrael & FelliHl978h : Gsi lKurtz et alJlIlQ^ : iui lSchneider et al.l ll2003h . 



Right ascension 


Right ascension 


Fig. 1 . Pointing for S 106 (left) and IRAS 23133+6050 (right). Overlaid on each GLIMPSE (iBeniamin et al.l2003HChurchwell et al.l 
l2009l) IRAC [8.0] image is the corrected SPIRE beam aperture (43" and 37" respectively). We adopt the convention of coloring 
points relevant to S 106 and IRAS 23133+6070red and blue respectively, beginning with the apertures in this Eigure. 


resolution radio/sub-mm observations (e.g.. [Richer et al.lll993h 
dispelled this notion though, and replaced it with the idea that the 
central regions of S 106 are dominated by a variety of clumps of 
denser material, which have in turn been s uggested to be de bris 
from the dissolution of a large dense torus (iLittle et al.lll995h . In 
subsequent years it has been suggested that the dark lane is the 
projected shadow of an edge on accretion disl0 around the cen¬ 


tral young star (e.g.. lNoel et al.l2()()5ll . howev er recent SOEIA re - 
sults seem to be inconsistent with this model (ISimon et al.l2012ll . 

Subsequently several authors have investigated the PDR 
environment of S 106 each using different datasets (e.g.. 


^ Other authors have detected the signature of such a disk without 
com menting on the issue of the projected shadow / dark band structure, 
e.g.. lMurak:awa et aT] ll201.3h . 
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van den Ancker et al.l l200(i 

ISchneider et alJ 20031 Noel et al.l 

20051). The first of these, by 

van den Ancker et al.l (120001). used 


IClegg et^ll996ll to find the global properties of the FIR emis¬ 
sion from S 106. They conclude, on the basis of very strong 
[Si ii] emission, that the PDR associated with S 106 must have 
densities of around 10^ cm~^ and UV radiation field strength 
of greater than 10^. Ivan den Ancker et alJ (l2000l) also conclude 
that this PDR mu st be 5-15" away from the central star. A sub¬ 
sequent study by ISchneider et all (l2003h showed that the PDR 
of S 106 is much more extended than previously thought, ex- 
tend ing up to several arcmi nutes from the central star. In addi¬ 
tion, ISchneider et al.l (l2003h used PDR diagnostics to show that 
a two phase PDR model was necessary to understand the line 
emission from S 106. They concluded that a very dense molec¬ 
ular component with an H 2 density greater than 3x10^ cm“^ 
was necessary, along with a less dense interclump phase. This 
finding agreed with the previous study of CO emission by the 
same authors which fo und evidence of a very dense component 
(ISchneider et al.l2002h . Parallel efforts to those using higher spa¬ 
tial resolution mid-IR inst ruments have yie lded further insights 
as to the central structures. iNoel et^(l2005l) studied S 106 using 
near-infrared spectro-imaging of a combination of hydrogen re¬ 
combination lines and molec ular hydrog e n emi ssion. From the 


_ 3ge n emi s! 

perspective of this work, the iNoel et aTI (l2005ll study came to 
two key results: 1) the observations of H 2 lines were found to 
be consistent with shocks driven by outfiows of the central star 
and, 2) the internal structure of the S 106 region based on their 
line measurements (presented in their Fig. 12). Shocks would be 
an expected part of any high-mass star forming region, however 
previous studies found litt le evidence for their presence (e.g., 
Ivan den Ancker et al.ll200(ih . 


3. Observations, Data Reduction and Anaiysis 

3.1. Observations 

A log of all Herschel observations used is presented in Table |2] 
The observations were performed as part of Herschel program 
OTl_epeeters_01. The objects were chosen to provide a sam¬ 
ple to detect the long wavelength PAH features. This paper con¬ 
cerns the PDR properties of those objects which were not in¬ 
cluded in other PDR surveys. 

The SPIRE observations each used the pointed spectroscopy 
instrument mode with no mapping to fill in the spatial gaps be¬ 
tween the pixels. The PACS observations used the spectroscopy 
‘Range Scan’ mode, which provides a complete spectrum from 
50-200 pm without any spectral gaps. These observations were 
combined to produce a complete spectrum for each object cover¬ 
ing from 50-650 pm for the area on the sky covered by the cen¬ 
tral SPIRE spaxel beam after correcting for constant beam size. 
The pointings for each object are shown in Eig.[T] along with the 
corrected SPIRE beam size for each source (see Section I3.2b . 
The spectral resolution varied as a function of wavelength for 
each instrument, with SPIRE providing resolutions between 
~400 and ~1300 and PACS giving resolutions between 1000 and 
5000. In practice this means that none of the detected lines were 
spectrally resolved. 

3.2. Data Reduction 

Spire The SPIRE data were reduced using the Hers chel 
Intera ctive Processing Environment (HIPE; version 11.0; lOtj 
I 2 OIOI) . Each of our observations was taken in the ‘Bright Source’ 


mode, which enables SPIRE to observe sources which would 
usually saturate the detectors. The reduction process used the 
“SPIRE Spectrometer Single Pointing User Reprocessing” HIPE 
script. This allows the user to calibrate bright source mode ob¬ 
servations and use the newest version of the SPIRE calibration 
data which at the time of reduction was SPIRE_CAL_11_0. 

Accurate calibration of our sources was achieved by using 
the HIPE semi-extended c orrection t ool w hich implements the 
methodology described bv lWu et'aP (l2013l) . The semi-extended 
correction tool calculates the degree to which the observed 
source differs from a point source using the break between the 
orders of the SPIRE detectors where the beam size changes. In 
order to use this tool it was necessary to reduce the observations 
using the ‘point source’ calibration rather than the ‘extended’ 
calibration. The semi-extended correction tool results in a spec¬ 
trum which has been corrected to a constant beam size, 43" 
for S 106 and 37" for IRAS 23133-H6050. This eliminates the 
fiux jump between the two orders of the SPIRE instrument and 
allows easier cross calibration with the PACS data. The semi- 
extended correction tool only corrects the fiuxes for the central 
pixel as the beams of the long and short wavelength detectors 
are centered with respect to each other. Eor each of the sur¬ 
rounding, off-axis, SPIRE pixels for both objects also contain 
semi-extended structure (e.g., Eig.[T]) and as such they cannot be 
reliably calibrated at this time due to this software limitation. 


PACS The PACS data were reduced using HIPE version 11.0 
and the PACS calibration set PACS_CAL_56_0. The observations 
of both S 106 and IRAS 23133-1-6050 were performed using the 
unchopped spectroscopy mode with corresponding background 
observations (see Table |2]i because both sources are extended 
beyond the maximum PACS chop throw. Exposure times and 
ObsIDs are given in Table|2] 

The PACS data are affected ^ several systematics, mainly 
order leakage and spectral ghostfl The regions of order leakage 
were trimmed as necessary. This principally affected the longer 
wavelength spectra in band R1 (105-210 yum), where it proved 
impossible to have a well calibrated spectrum beyond around 
190 /rm. Regions affected by these problems were trimmed from 
the final data, resulting in gaps in the spectral coverage. 

It was found that the continuum was offset slightly in each 
PACS spectral order, likely a refiection of the varying beam size, 
and as such the orders were scaled such that they agreed with 
the longest wavelength order. It was found that the offsets be¬ 
tween the orders were additive, in that the line fiuxes for lines 
in overlapping regions matched before scaling. For S 106, the 
blue orders agreed with the final red order (B2A, B2B, RIB. 
Hence the only scaling necessary was for R1 A, which had to be 
increased by the addition of 850 Jy to agree with the continua of 
the other orders. For the IRAS 23133-1-6050 measurements, the 
first PACS spectral order was found to have a higher background 
than expected resulting in a strong offset of 950 Jy with respect 
to the continuum of the other bands. While these values sound 
large, they represent factors of only around 20%. 

The final PACS spectra were created by coadding the pixels 
to match size of the beam produced by the SPIRE semi extended 
source correction. This was accomplished by calculating the rel¬ 
ative contributions of each pixel to this beam and summing all 
of the pixels in the correct proportions. The final fluxes were 
achieved by scaling this coadded PACS spectrum such that the 
continuum slopes and intensities matched. This scaling factor 

see the PACS data processing guide: 
http://herschel.esac.esa.int/Data_Processing.shtml 
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Table 2. Observation Log 


Object 

Date 

Instrument 

Mode 

Obs. ID 

Duration (s) 

A Range (pm) 

Comment 

S 106 

2011-05-25 

SPIRE 

Spectro Point 

1342221683 

528 

190-650 


S 106 

2011-11-20 

PACS 

Range Scan 

1342232575 

670 

50-74 & 100-145 


S 106 

2011-11-20 

PACS 

Range Scan 

1342232576 

339 

50-74 & 100-145 

background 

S 106 

2011-11-20 

PACS 

Range Scan 

1342232577 

669 

70-105 & 140-210 

background 

S 106 

2011-11-20 

PACS 

Range Scan 

1342232578 

1330 

70-105 & 140-210 


IRAS 23133 

2011-08-15 

SPIRE 

Spectro Point 

1342227453 

1068 

190-650 


IRAS 23133 

2012-01-14 

PACS 

Range Scan 

1342237480 

670 

50-74 & 100-145 


IRAS 23133 

2012-01-14 

PACS 

Range Scan 

1342237481 

339 

50-74 & 100-145 

background 

IRAS 23133 

2012-01-14 

PACS 

Range Scan 

1342237482 

669 

70-105 & 140-210 

background 

IRAS 23133 

2012-01-14 

PACS 

Range Scan 

1342237483 

1330 

70-105 & 140-210 



was found to be multiplicative in order for the fluxes of the '^CO 
lines to be consistent across the PACS/SPIRE gap. For S 106 this 
factor was relatively small (1.5), while for IRAS 23133+6050 it 
was much greater (4). 


3.3. Analysis 

3.3.1. Line Fluxes 

In each of the coadded final spectra, the spectral lines were an¬ 
alyzed in the following ways. The line fluxes were extracted via 
fitting Gaussians and local continua to the data. It was found that 
better fits resulted from fitting a linear local continuum to regions 
of the spectrum on either side of the line, and then fixing the lin¬ 
ear parameters fitted for these regions when fitting the full spec¬ 
trum of a Gaussian line profile plus linear continuum. The signal 
to noise estimates for each line measurement were calculated by 
using the root mean square (rms) noise in the continuum spec¬ 
tral windows after the subtraction of the linear fit. This rms value 
can then be combined with the integrated line flux, width and the 
instrumental resolution to yield the final signal to noise ratio for 
the measured line. There is only one significant blended pair of 
lines present in the data, the 371.65 fj.m CO 7-6 and 370.42 //m 
[C i] 2-1 lines. Initially we performed a fit using two gaussians to 
the complex. We have checked this fit by calculating the flux of 
the 7-6 CO line using the neighboring CO lines and subtracting 
this flux from the total flux of the blended complex. The residual 
of this process yielded the same [C i] flux as the initial fitting 
procedure. The measured line fluxes are given in Table [3] 

For S 106 we can compare the measured PDR cool¬ 
ing line fluxes here to those found by other observations in 
the past. The most direct compari s on ca n be drawn to those 
presented by Ivan den Ancker et al.l (l2000t) and ISchneider et al.l 
(l2003l) which are based on observations by IS O-LWS, KAO- 
FIFI and KOSMA-SM ART (IClegg et akllio^ iPoglitsch et all 
Il99lt iGraf et ani2002h . A direct comparison between the mea¬ 
sured surface brightnesses is given in Table|4] Some care should 
be taken in comparing the presented quantities as the beam 
sizes differ by large factors from the Herschel measurements 
presented here. For the dominant cooling line, [O i] 63 yum, it 
is easy to see why the ISO-LWS surface brightness measure¬ 
ments would be much smaller as they were measured using a 
larger beam by a factor of around five, and hence the over¬ 
all surface brightness is lower as the flux of S 106 is peaked 
strongly in the central regio ns near the exciting source IRS4. 
Ivan den Ancker et akl (l2000l) measured surface brightnesses for 
a separate set of ISO-LWS observations which were also cen¬ 


tered on S 106 and found that the [O i] 6 3 um line had a s urface 
brightness between the values found by ISchneider et al.l (l2003h 
and our measurements. Evidently different pointings around the 
S 106 region can result in very different results, even consider¬ 
ing the fact that the comparatively large LWS beam encompasses 
the majority of strong PDR cooling line emission from the S 106 
complex (e.g. see the Sofla/G REAT map of [C ii] 158 fim pre¬ 
sented by Simon et"^ (1201^ 1. 

Overall, our observations of the [O i] line agree within ~30% 
with the LWS measurements if we assume that the S 106 flux is 
concentrated at the center of the beam and compensate for the 
different beam sizes. The same process yields similar results for 
the other LWS cooling line measurements - for example for the 
[C It] 158 jum line ou r results agree with th e LWS observations 
almost exactly for the lyanjen_Ancker et aljj [2000l) pointing and 
to within 30% for the ISchneider et al~ (l2003l) po i nting, assum¬ 
ing that the beam size quoted bv ISchneider et al.l (l2003l) is cor¬ 
rect (see footnote to Table Ell. Overall then, we find agreement 
within the quoted absolute flux uncertainties of 30% for both 
ISO-LWS and PACS ILWS: ISwinvard et al.lll 9961: PACS; PACS 
Spectrometer Calibration Documen{3)- At first glance, the KAO 
surfac e brightness measurement of [O i] (from Schneider et al.l 
l2003l) that we include in Table El seems incompatible with our 
measurement, as it was for a beam smaller than the SPIRE 
beam by a factor of three and fo und a lower surface brightness. 
However. lSchneider et al.l (l2003l) discuss this problem as the dis¬ 
crepancy between KAO/Herschel i s almost the same as bet ween 
the KAO and LWS measurements. I^hneider et al.l (l2003l) con¬ 
clude that the KAO and LWS measurements can be reconciled 
by taking into account the spatially varying emission intensities 
- specifically that if one assumed a gaussian profile for the [O i] 
emission the LWS beam size, the LWS flux could be recovered 
from the KAO measurement. 


3.3.2. Total FIR Flux 

The total flux of the far-infrared dust continuum was determined 
by fitting a modified blackbody to each spectrum. In order to 
better constrain the peak of this spectrum, we searched for addi¬ 
tional observations of our sources in the 20-60 pm regime such 
that we might acquire extra data points via aperture photome¬ 
try using the SPIRE beam. However, realistic measurements of 
the fluxes for these sources were not available as they are very 
bright and saturated in both MIPS and WISE data. Figure l^in- 


^ http://herschel.esac.esa.int/twiki/pub/Public/PacsCalibrationWeb/ 
PacsSpectroscopyPerformanceAndCalibration_v2_4.pdf 


5 


























































D. J. Stock et al.: Herschel spectroscopy of Galactic PDRs 



10 100 1000 

Wavelength, 


Fig. 2. Full PACS + SPIRE spectra for IRAS 23133+6050 (top, blue) and S 106 (bottom, red, multiplied by 0.01 for clarity). In 
addition we include the IRAS photometric points for each sourc e (triangles), and for IRA S 23133+6050 also the WISE photometry 
dWright et al.ll2010h as squares and the SCUBA observations of lThomnson et al] (l2006h as diamonds. A modified blackbody fit for 
each spectrum is shown as a dashed line. 


eludes IRAS photometric points, however these are integrated 
totals for the entire source and not scaled to the Herschel beam. 
In addition we investigated using MSX photometry, however the 
MSX fluxes agree well with the IRAS points. We therefore in¬ 
vestigated using spectral data from ISO-SWS and LWS to help 
constrain the shorter wavelength regime. The ISO spectrome¬ 
ters used very different beam sizes though, with the LWS aper¬ 
ture being of ~ 120" diameter while the SWS aperture was 20" 
X 33"at its longest wavelengths. This corresponds to very dif¬ 
ferent areas of each object dominating the beam. Eor SWS the 
bright central parts of the nebulae All the beam, while for LWS 
the colder outer areas All the beam. For S 106, the LWS spectra 
agree well with the Anal PACS spectrum, albeit with an ofl'set 
of around 50%. This is likely because there is extended bright 
structure in the LWS beam and as such it gives a similar spec¬ 
tral slope. Lor IRAS 23133+6050, the picture is less clear, as 
the PACS spectrum has a much steeper gradient than the LWS 
spectrum. IRAS 23133+6050 is much more compact than S 106, 
with all of the hot material being conflned to the center. The shal¬ 
lower gradient of the LWS continuum then may reflect the extra 
cold material in the beam. 

Lor each object the peak of the SED falls just before the 
PACS wavelength range, therefore in order to measure the to¬ 
tal infrared flux we fit a modified blackbody to the spectrum and 
integrated all of the flux of the fitted blackbody between 10 and 
1000 //m. The aperture mismatches with the existing ISO-SWS 
and -LWS data and our observations proved impossible to rec¬ 
oncile and as such we have no strong constraints on the blue 


end of the spectrum. As a test of this uncertainty, we performed 
an additional fit in which the FIR flux was fit with two black- 
bodies. The single blackbody fit discussed previously was held 
constant, while an extra component was added to fit the IRAS 
points shortward of the Herschel data. The FIR flux found in this 
way was around 25% greater than that measured using only one 
blackbody. This figure is likely to be an overestimate of the true 
FIR flux, as the flux of the 12 pm IRAS Alter is dominated by 
PAH emission. For repeatability and simplicity, we adopted the 
single blackbody measurement with an uncertainty of 25% on 
the derived total FIR fluxes. The total FIR fluxes for each object 
derived in this manner are listed in Table[3] 


4. PDR Diagnostics 

It is clear from inspection of Table [2 that the [O i] 63 pm line 
dominates cooling for both PDRs, dominating both of the other 
important cooling lines ([O i] 146 pm and [C ii] 157 pm) by fac¬ 
tors of between five and eight. This is typical of dense PDRs. 
In addition CO rotational lines extend to very high J levels also 
indicating warm dense PDRs. On the other hand, the [C ii] and 
[C i] emission can also arise in low density gas illuminated by 
weak radiation fields. There is likely a mix of densities, temper¬ 
atures, and Ailing factors that combine to produce the observed 
emission. We next use PDR line diagnostics to investigate the 
various phases. 
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Table 3. PDR Line and Continuum Flux Measurements 


Species 

Transition" 

Wavelength 

[pm] 

Flux 

[10^^ erg s“* sr“’] 

IRAS 23133-H6050 S 106 

“'^CO 

4- 

3 

650.25 

5.67 ± 0.43 

6.96 ±0.16 

‘^CO 

5- 

4 

520.23 

8.11 ±0.18 

10.6 ±0.1 

'^CO 

6- 

5 

433.56 

11.6 ± 0.1 

14.3 ±0.1 

‘^CO 

7- 

6 

371.65 

16.4 ± 0.2 

20.1 ±0.3 

‘^CO 

8- 

7 

325.23 

20.6 ± 0.6 

26.6 ± 0.5 

‘^CO 

9- 

8 

289.12 

28.1 ±0.3 

29.6 ± 0.5 

‘^CO 

10- 

-9 

260.24 

35.6 ± 0.5 

39.2 ± 0.5 

‘^CO 

11 - 

10 

236.61 

37.7 ± 0.6 

40.4 ± 1.0 

‘^CO 

12- 

11 

216.93 

45.0 ± 1.0 

50.5 ± 1.7 

‘^CO 

13- 

12 

200.27 

45.7 ± 2.5 

54.1 ± 3.4 

‘^CO 

14- 

13 

186.00 

47.2 ± 0.8 

57.3 ± 0.5 

‘^CO 

15- 

14 

173.63 

42.4 ± 0.7 

49.5 ± 0.4 

‘^CO 

16- 

15 

162.81 

35.8 ± 1.4 

44.5 ± 2.4 

‘^CO 

17- 

16 

153.27 

28.9 ± 1.0 

35.9 ± 0.9 

‘^co 

18- 

17 

144.78 

20.6 ± 2.5 

24.1 ± 1.7 

‘^co 

19- 

18 

137.20 

17.7 ± 1.6 

19.1 ±0.4 

‘^co 

20- 

19 

130.37 

12.9 ± 1.0 

13.6 ±0.4 

‘^co 

21 - 

20 

124.19 

6.43 ± 2.00 

11.7 ±0.5 

‘^co 

22- 

21 

118.58 

4.11 ± 1.53 

10.3 ±0.8 

‘^co 

23- 

22 

113.46 

4.34 ± 2.97 

8.74 ± 0.89 

‘^co 

5- 

4 

544.16 

1.83 ±0.13 

1.81 ±0.21 

‘^co 

6- 

5 

453.50 

2.49 ±0.13 

2.91 ± 0.06 

‘^co 

7- 

6 

388.74 

2.14 ±0.18 

2.47 ±0.14 

‘^co 

8- 

7 

340.18 

3.66 ± 0.42 

5.52 ± 0.49 

‘^co 

9- 

8 

302.41 

5.53 ± 1.57 

5.37 ± 1.38 

‘^co 

10- 

-9 

272.20 

4.59 ± 0.45 

6.77 ± 0.91 

‘^co 

11 - 

10 

247.49 

5.29 ± 0.63 

6.05 ± 0.83 

‘^co 

12- 

11 

226.90 

3.47 ± 0.88 

4.45 ±1.10 

[C I] 

^Pi- 

'Po 

609.14 

1.36 ±0.09 

0.86 ±0.14 

[C I] 

^P2- 

^Pi 

369.87 

3.09 ± 0.23 

2.72 ± 0.27 

[Cii] 

- 


157.74 

604.7 ± 3.5 

391.1 ±2.0 

[01] 

^Pl - 

^Po 

145.54 

439.6 ±2.1 

447.9 ± 1.7 

[Nil] 

3P2- 

'Pi 

121.91 

34.3 ±2.1 

34.7 ± 1.3 

[0 III] 

^Pl- 

'Po 

88.36 

34.2 ± 4.6 

233.1 ±2.3 

[01] 

^P2- 

'Pi 

63.18 

2934 ± 10 

2987 ± 7 

EFIR 

- 


10 - 1000 

1.86 X 10« 

1.69 X lO*^ 


Upper and lower J values for *^CO and *^CO. 


4. 1. FIR Continuum and Stellar Spectral Types 

The FIR continuum can be used to directly probe the UV radi¬ 
ation field impinging on the PDR system. The total energy ab¬ 
sorbed by the dust in the syste m must be the same as that re¬ 
emitted. Previous authors (e.g., iMeixner et al.lll9^ have used 
this relationship to calculate a va lue for Gn. This calc ulation 
is summarized in Appendix B of IMeixner et al.l (Il992h where 
it is shown that the derived Go values weakly depend on the 
emitting structure. For a spherically symmetric source (e.g., 
IRAS23133-I-6050), it can be assumed that the absorbing and 
radiating surface areas are equal, and hence the relationship can 
be summarized as: Go = 2FIR / 1 Habing field. For more com¬ 
plex sources the relationship between the r adiating and absorb- 
ing surface areas is different, for example. IMeixner et al.l[l9^ 
show that for a disk system that the derived Go will be different 
by a factor of two. S 106 would seem to be a mixture of these 
two paradigms, with a small disk around the central star being 


further encased within a molecular cloud/PDR, therefore there 
will exist a range of values of Go spanning these geometries. For 
each object the resulting values of Go assuming a spherically 
symmetric geometry are very similar at around 1.3x10^. For the 
disk region of S 106 though. Go could be up to 2.6x10^. 

An additional measure of Go can be calculated using the dis¬ 
tance from the ionizing source to the PD R and the stellar spectral 
type of the exciting star. iTielensI (1201 Oh gives the relationship be¬ 
tween Go and the stellar luminosity as: 

Go - 625 X ^ (1) 

47rfl“ 

where L, is the stellar luminosity (in ergs s^'),;^^ is the frac¬ 
tion of photon energy emerging from the H ii region between 6 
and 13.6 eV, and d is the distance between the star and the ioniza¬ 
tion front in cm. Both IRAS 23133-1-6050 and S 106 are powered 
by late 08-9 type stars, which have far UV luminosities around 
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Table 4. Prior Measurements of PDR cooling lines in S 106 


Line 

Wavelength 

Herschel Brightness 

Herschel Beam 

Prior Measurement 

Beam Instrument 


[pm] 

[10^^ erg s“' cm“^ sr“*] 

[10-® sr] 

[10^^ erg s“' cm“^ sr“*] 

[10-® sr] 


vandenAncker_elal^(200^: 


[O I] 

63 

2987 ± 7 

3.41 

1125 ± 114 

16.3“ 

ISO-LWS 

[0 I] 

63 

2987 ± 7 

3.41 

1125 ± 114 

3.4U 

ISO-LWS 

[0 I] 

146 

447.9 ± 1.7 

3.41 

288.8 ± 29.0 

8.8“ 

ISO-LWS 

[Cii] 

158 

391.1 ± 2.0 

3.41 

157.2 ± 15.9 

11.6“ 

ISO-LWS 


Schneider^el^l^(200^: 


[O I] 

63 

2987 ± 7 

3.41 

639 + 7 

12.4 

ISO-LWS 

[O I] 

63 

2987 ± 7 

3.41 

2100 

1.3 

KAO-FIFI 

[O I] 

146 

447.9 ± 1.7 

3.41 

74 + 3 

9.0 

ISO-LWS 

[Cii] 

158 

391.1 ± 2.0 

3.41 

80 

8.6 

KAO-FIFI 

[Cii] 

158 

391.1 ± 2.0 

3.41 

107 + 3 

8.5 

ISO-LWS 

[C I] 

609 

0.86 + 0.14 

3.41 

0.217 

7.1 

KOSMA 

[C I] 

370 

2.72 ± 0.27 

3.41 

2.64 

7.1 

KOSMA 


°: ISchneider et all ( 120031) quote smaller values for the area of the LWS beam at the same wavelengths than ivan den Ancker et all ( |200(]1) . which is 
likely due to the ISO-LWS beam profiles being significantly different in practi ce than expected, le ading to changes to the quoted beam sizes in the 
LWS documentation o ver time. For example see the discussion presented bv iLiseau et al.l l l2006l . Section 3.1). As such, the values presented by 
ISchneider et al.l (1200 3ll are likely to be closer to the true LWS beam sizes. 


3.25 X lO'* times that of the suifl. For IRAS 23133+6050, the 
ionized region has a radius of 4.5", which at a distance of 5.5 
kpc translates to 0.1 pc. This implies a Go value of around 5 
X 10"^. Similarly for S 106, the radius of the ionization front is 
aroun d 0.1 pc (based o n a distance of 1.5 kpc and a radius of 
12.5": iNoel et al.|[2005h . which ir nplies a Go of about 8 x 10"^ 
when using the same assumptions. iNoel et al.l (l2005h found that 
the H II region of S 106 was approximately cylindrical, and as 
such the Go we calculate using their radius is only appropriate 
for the parts of the ionization front closest to the exciting star 
and therefore is actually an upper limit. In good agreement with 
our analysis, the Habing field for S 106 determined from the 
Herschel PACS 70 and 160 micron fluxes is 3 x 10"^ Go(value at 
a distance of 0.1 pc from S106 IR in a 12” beam, Schneider et 
al., in prep). 

4.2. Diagnostic Diagrams 

The measured line fluxes given in Table [3 can be combined in 
various ways to produce ratios which are sensitive to the physical 
conditions in the PDR. The variation of pairs of these ratios with 
respect to the fundamental PDR parameters Go and iflare shown 
in Figs.[3]and|4l To c onstruct these diagra ms we used the results 
of the PDR models of lWolflre et al.l (1201 Ol) and iHollenbach et alJ 
(l2012h . There are a wide variety of different PDR models avail¬ 
able in the literature, and they have been extensive ly compared 
and va lidated against each other, for example see iRollig et aH 
(|2^, in which the model we have chosen to employ here is 
referred to as the HTBKW model. These models have all been 
shown to reasonably agree for most PDR diagnosti cs and all re¬ 
produ ce the same trends when run on idential data (iRollig et alJ 
1200 ^ . 

In each case we consider a range of Go values between 10^ 
and 10^ Habings and a range of densities between 10^ and 10^ 


® iParravano et all ( 120031) provide the relationship between stellar 
mass and the FUV flux, so the quoted figure, which represents L,x in 
Equation[T] can be arrived at by assuming a stellar mass of around ~21 
Mq for an 08-9 star. 


cm~^. In Table |5] we give an overview of the physical condi¬ 
tions indicated by the various diagnostics discussed in addition 
to those implied by the total FIR flux and the stellar spectral 
types of both sources. 

4.2.1. [C ll] AND [Oi] 

The primary diagnostics usually employed for PDRs involve ra¬ 
tios between FIR fine-structure cooling lines. Usually the ratios 
involving [C ii] 158 fim and the two [O i] lines are employed 
reflecting that these lines are the best understood of those in¬ 
cluded in PDR models. Typically such ratios do not correspond 
to unique PDR parameters, as their dependence of the line ra¬ 
tios on Go and n is complex. In order to produce a single value, 
other diagnostics must be employed such that intersections at 
specific Go and n values are found. The simplest combination of 
diagnostics which yields such physical conditions occurs when 
considering [C ii] 158 /rm / [O i] 63 gm and ([O i] 63 ;um + [C ii] 
158 fim) / SFIR diagnostic diagram (Fig. [3 top). 

This diagram implies that both PDRs have similar physical 
conditions, with log Go around 4 and log n around 3.75, well 
within typical PDR conditions. These values are compatible with 
those which were discussed in Section|3 For IRAS 23133+6050 
the density of ~ 10^ agrees nicely with that expected given the 
density of the H ii region, while for S 106 this density likely rep¬ 
resents the lower density material which make s up the majority 
of the PDR volume ('e.g.. lSchneider et^l2003h . 

The ratio of the [O i] lines measured for both sources is of the 
order of [O i] 145 um/ [O i1 63 um ^ 0.15 , which agrees with the 
measurements of ISchneider et alJ (l2003h who found 0.17 using 
ISO-LWS observations. This clearly ind icates that one or both 
of the [O i] lines is optically thick (e.g.- lTielens & Hollenbachl 
Il985bl : Fig. 2): if both lines were optically thin their ratio would 
be less than 0.1. More recently it has been shown that in many 
cases the [O i] line ratio indicates optically thick conditions, 
though other explanations for such ratios exist. It was concluded 

’ Henceforth we use n to refer to the density of hydrogen nuclei such 
that for a purely molecular gas n(H 2 ) = 0.5 x n. 
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Table 5. Results of Simple PDR Diagnostics" 


Diagnostic 

Figure 

log Go 

log n 

log Go 

log n 



IRAS 23133+6050 

S 106 

General: 






Spectral Types 


4.6 


14.9 


SFIR 


5.1 


5.1-5.4 


Low density 






[C ii] 158 jim __ [C ii] 158 jim + [O i] 63 jim 
[O i] 63 um ’■" zFIR 

m 

4.0 

3.5 

4.2 

3.8 

[C II] 158 pm [C II] 158 pm + [0 i] 145 pm 
[0 i] 145 pm ’’" eFIR 

m 

4.0 

4.1 

4.2 

4.5 

[Ci] 370pm [Cii] 158pm+[Oi] 63 pm 
[C i] 609 pm iFIR 

m 

3.2 

2.5 

3.6 

3.2 

High density 






■-CO (J=7-6) [C ii] 158 pm + [0 i] 63 pm 
12CO (J=4-3) eFIR 

a 

4.9 

4.7 

4.9 

4.7 

“CO (1=15-14) 

“CO (1=7-6) 

a 

5.2 

5.6 

5.2 

5.6 

[C ii] 158 pm + [O i] 63 pm 






“CO(l=15^if 

“CO (1=4-3) 

a 

5.2 

5.4 

5.2 

5.4 

[C ii] 158 pm + [0 i] 63 pm 






eFIR 






Questionable: 






[C i] 609 pm [C II] 158 pm + [O i] 63 pm 
“CO (1=4-3) eFIR 

m 

5.2 

5.4 

5.2 

5.8 


With Go expressed as a multiple of the Habing Field and n in cm 


by iLiseau et alJ (l2006h that self absorption of the [O i] 63 ;um 
line along the line of sight with modest optical depths of around 
1-2 could be resTOnsible for high values of the [O i] 145 /tm / 
[O i] 63 fim raticfl Those results refer mainly to an ISO-LWS 
sample of dark clouds, for which the density required for both 
[O i] lines to be optically thick was considered very unlikely. 
For IRAS 23133+6050 and S 106 though, such densities may 
be plausible - however [O i] 63 //m self absorption would likely 
also play a role in shaping the flux ratio, rendering it difficult to 
disentangle the two. 

4.2.2. [C I] 

In addition, we also investigate the [C i] hne-structure lines by 
considering the ratio of the fluxes of the 370 /rm line over the 
609 jum line (bottom panel of Fig. [3]l. This ratio indicates very 
different physical properties to those mentioned previously, with 
the ratio of the [C i] lines indicating low density, weak radiation 
field PDR conditions. It should be noted though that it has been 
shown that the [C i] intensities given by models are generally un¬ 
certai n (maximal variation of around a factor of ten: lR611ig et alJ 
l2007h and therefore the results for this line ratio are highly model 
dependent. 


* ILiseau et all ( 120061) actually consider the inverse of this quan¬ 
tity and as such their higher/lower than adjectives are switched, we 
adopt the same sens e as is provided by the PDRT diagnostics of 
iKaiifman et Zl (l2006h : IPound & Wolfirel d2008l) . 


4.2.3. 

In addition to the hne-structure cooling line diagnostics we can 
also employ the pro perties of the molecul ar lines in the spec¬ 
tra in the same wav. iKaufman et ^ (1 19991) provide the various 
diagrams for the different '^CO transitions, here we investigate 
ratios of t hree: 4-3, 7-6 and 1 5-14 , but using the updated P DR 
models of lWolflre et al.l (1201 Ol) and Hollenbach et al.l (l2012l) . 

The ratio of the two mid-J lines: *^CO(J=7-6)/'^CO(J=4- 
3) is shown in Fig. |4] (top), and indicates a much stronger UV 
field (log Go = 4.9) and denser environment (log n - 4.8) for 
both sources than the conditions indicated by the fine-structure 
lines. For the highestj'^C01ine,J=15-14, even denser and more 
strongly illuminated conditions are found, with densities an or¬ 
der of magnitude higher and radiation fields around double that 
found from the lower excitation lines. Diagnostic diagrams for 
the (1=15-14) line, as compared to the lower excitation CO 
lines, are shown in the lower two panels of Fig. 01 

In addition to the pure ratios, we also show the ratio of 
the [C i] hne-structure line with the '^CO(J=4-3) line in Fig.|5] 
This diagnostic clearly shows a similar pattern to that of the other 
'^CO diagnostic ratios, i.e. high log Go (~ 5) and log n (~ 5). 

4.3. CO Ladder fluxes 

In the previous section we exploited the simple relationships be¬ 
tween different CO transitions resulting from PDR models to 
And the physical conditions. It is also possible to derive such 
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Fig. 3. PDR diagnostic diagrams fo r observed fi n e-stru cture 
lines - based on the PDR models of IWolfire et al.l (1201 Oh and 
iHollenbach et alJ (l2012h . Top: [C ii] 158 pm / [O i] 63 pm ver¬ 
sus ([O i] 63 pm + [C It] 158 pm)l ZFIR. Middle: [C ii] 158 / 
[O i] 145 versus ([O i] 145 pm + [C ii] 158 pm)l SFIR. Bottom: 
[C i] 370 / [C i] 609 versus ([O i] 145 pm + [C ii] 158 pm)/ 
SFIR with the density grid extended down to log n = 2. In each 
figure the blue and red points represent the values found for 
IRAS 23133H-6050 and S 106 respectively. 


Fig. 4. PDR diagnostic diagrams for ratios of '^CO(J=15-14), 
'2 cO(J= 7-6) and '2cO(J=4-3) versus ([O i] 63 pm + [C ii] 158 
pm)/ SFIR. In each figure the blue and red points represent the 
values found for IRAS 23I33H-6050 and S 106 respectively. The 
data points indicate a much denser environment with a stronger 
UV field than that seen for the fine-structure line diagnostics. 
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Fig. 5. [C i] 609 / '^CO (J=4-3) versus ([O i] 63 //m + [C ii] 

158 fim)! SFIR PDR diagnostic diagram. The blue and red points 
represent the values found for IRAS 23133+6050 and S 106 re¬ 
spectively. 



E, K 


Fig. 6. '^CO (squares) and (triangles) ladder fluxes for 

IRAS 23133+6050 and S 106 in terms of the excitation energy 
of that level. 


conditions directly from the fluxes of the whole set of molecular 
lines. In Fig. we show the and '^CO line fluxes against 
the energy of the upper level Ey. The detections of the '^CO 
lines for both IRAS 23133+6050 and S 106 are measured across 
the entire SPIRE/PACS range, that is, from J„ = 4 to 23, while 
for '^CO the range is smaller, from J„ = 5 to 12, owing to their 
weaker fluxes. The '^CO fluxes shown in Pig.|6](a) are remark¬ 
ably similar for both objects, peaking at around Ej = 550 K, or 
J„ = 13 or 14. The only difference between IRAS 23133+6050 



Fig. 7. Combination of RAD EX models scaled to fit the S 106 
'^CO ladder. The low density component with log n - 4, T = 
300 K is represented by the green line, while the dominant con¬ 
tribution arises from the dense component with log n - 5.6 and 
a temperature of 400 K which is represented by the purple line. 
The aggregate model (black line) matches the observed '^CO 
ladder from Jy = 4 to 20. 


and S 106 comes from the higher J lines, where for S 106 the 
four highest J lines show additional excitation, whilst those for 
IRAS 23133+6050 continue to decrease. The possible implica¬ 
tions of this will be discussed in the following section. The '^CO 
line fluxes seem to peak at around E^ = ~300 K, or J„ = 10, how¬ 
ever the data are noisy and it is not clear that this apparent peak 
is real. 

The peak of the '^CO fluxes is itself an indicator of 
the physical conditions of both sources. Using the density of 
the high density component predicted by the simple diagnos- 
tics (n ~ 10^-^ cm~^) and the radiative transfer code RADEX 
(Ivan der Tak et al.ll2007h we calculated the fluxes of *^CO lines 
for a variety of temperatures. The required temperature to fit the 
peak is then around 400 K. In Eig|2]we show this model, com¬ 
bined with a model representing the low density component at 
the same temperature, which have been scaled to match the over¬ 
all flux of the whole CO ladder. However, the combination of 
density and temperature which fits the peak is not unique, and 
higher or lower density models can fit the peak almost equally 
well with suitably reduced or increased temperatures. Our choice 
of n = 10^ ® cm“^ is driven by the simple diagnostics discussed in 
previous subsections. In contrast, the highest and lowest ‘plau¬ 
sible’ solutions considered have densities and temperatures of 
10’ or 5x10^^ cm“^ and 200 or 1500 K respectively. Thus RADEX 
models lead us to conclude that the most plausible explanation 
for the mid-J (Jy/ = 6-20) *^CO ladder is a (very) high density 
PDR with a temperature of around 200^00 K. We attempted 
to fit the high-J S 106 lines using RADEX, however it was not 
possible to And a good fit to those points without severely com¬ 
promising the quality of the fit to the other points. 
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E, K 


Fig. 8. '^CO (diamonds) and '^CO (triangles) rotation diagrams 
for IRAS 23133+6050 and S 106. Two ultra-compact H ii re¬ 
gions from the iKarska et al.l (12014ll sample which display the 
warm CO component are also shown as orange, purple lines. 
Simple fits to the lower J CO lines are shown to demonstrate the 
deviation of the higher J CO levels from this trend for S 106 and 
the Karska sources. A fit to the highest J CO lines for S 106 is 
also included. The Karska sources have been offset in y for clar¬ 
ity. The first '^CO and '^CO points for both of our sources are 
4-3 and 5-4 respectively. 


4.3.1. CO Rotation (diagrams 

To further illustrate the differences between our sources and 
those found in the literature we create rotatio n diagrams from 
the *^C O and '^CO fluxes (Fig. [8]l. Following Llusttanont et al.l 
(I 2 OOOI) we construct rotation diagrams showing the relationship 
between the line fluxes Fji, the temperature and the column 
density No as follows: 


EjIkT, = In 


No \ 

‘ 5 ( 7 ’,)^-/ 


_47tF^\ 

AjihvjigjJ 


( 2 ) 


where gj refers to the statistical weight of the upper level 
j, Vji is the frequency of the transition, <h(7’x) is the partition 
function, Aj, is the oscillator strength f or the particular transi- 
tion and d is the distance to the source. Llusttanont et alJ (l200(ll) 
refer to the second term as Yj and we follow this convention. 
The temperature and column densities can then be found in the 
limit of optically thin emission in local thermodynamic equilib¬ 
rium (LTE) by fitting a straight line to the Ej, Yj coordinates. 
There are many potential problems with this method, mainly 
involving the aforementio ned assumptions (see, for example, 
[Goldsmith & Lang^ll999h . In light of that, it is not surprising 
that both IRAS 23133+6050 and S 106 do not follow exactly 
straight lines. Nevertheless, since it is a common procedure, we 
will analyze the observations using this technique and compare 
to other methods and models. 


Table 6. Properties of fits to Rotation Diagram Components 


Source 

Species 

Ju’s 

T[K] 

log(N) 

IRAS 23133+6050 

‘^CO 

all 

188 

49.33 

IRAS 23133+6050 

‘^CO 

all 

94 

48.84 

S 106 

‘^CO 

4-20 

190 

48.91 

S 106 

‘^CO 

20-23 

327 

48.00 

S 106 

‘^CO 

all 

96 

48.40 


The curvature exhibited at low Ej is expected and usually in¬ 
terpreted as requiring a second, warm, temperature component 
(c.f. Eigure |7]i. These levels trace the warmish material on the 
surface of the (putative) molecular cloud, the colder molecular 
cloud material would only manifest itself in unobserved lower J 
transitions (J j 3). However recent work has shown that this cur¬ 
vature can be explained in the optically thin limit by deviations 
from LTE at den sities below 10^ cm“^ and constant temperatures 
(lNeufeldll2012h . This effect arises because, in the low density 
limit, the population is dominated by the low lying transitions, 
with collisional excitation and spontaneous decay balanced for 
the higher level transitions. 

Erom these best fits to the *^CO lines of S 106 and 
IRAS 23133+6050 we find excitation temperatures of around 
190 K and around 10"^® molecules along the line of sight (sum¬ 
marized in Table|6]l. We have also included a linear fit to the high 
excitation component of the S 106 *^CO spectrum - this fit yields 
a much higher excitation temperature of around 330 K with a 
lower density. The derived excitation temperatures for '^CO are 
systematically lower than those found by the RADEX models dis¬ 
cussed in the previous subsection. This discrepancy is likely due 
to the assumptions about optical depth which underpin the use 
of rotation diagrams - it seems unlikely that many of the low 
J CO lines are optically thin, for exampl^l The RADEX models 
make no such assumptions and are therefore a somewhat more 
trustworthy indicator of the temperature of the various compo¬ 
nents, on the other hand the RADEX models assume either an ar¬ 
bitrary gas temperature distribution or constant temperature, and 
constant molecular abundances. However, we could not find a 
plausible solution for the very high excitation S 106 points with 
RADEX models, whereas the rotation diagrams tell us that the 
temperature of that component is higher than that of the medium 
density component by at least 100 K. However, even this argu¬ 
ment relies on the assumption that our sources are in LTE, which 
is not altogether clear from inspection of Figure [8] As such the 
parameters of the rotation diagram fits should be regarded with 
some caution, we will discuss how well these results compare to 
the results of the more detailed models in Section lOl 

The linear component seen in the four highest J *^CO lines 
of S 106 betrays the possible existence of an additional phys¬ 
ical component. There appear to be two possible explanations 
for this phenomenon: an even higher density PDR; or shocks 
in the S 106 environment. Such additional excitation in the up¬ 
per CO levels has been observed commonly in a related class 
of objects: protostars. We have therefore included in the rota- 

® A simple estimate of the *^CO optical depth can he made from the 
ratio, R, of the *^CO / '^CO line intensities. The optical (lepth, j n = 
X/R where X is the isotope a bundance ratio feurton et alll2013l) . For 
plausible values of X (e.g. see lMilam et alJl200.5h . the ‘^CO line is op¬ 
tically thick to the highest J level for which both *^CO and *^CO lines 
are available. 
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tion diagram several prominent examples from a recent survey 
of hig h mass star forming regions using PACS by iKarska et al.l 
(l20ll . The two objects (G34.26+0.15 and NGC 7538) were 
chosen for comparison on the grounds that they are ultra¬ 
compact H II regions and therefore closest in structure to S 106 
and IRAS 23133-H6050. In fact, NGC 7538 is actually a bipolar 
H II region similar to S 106. It is clear that these sources display 
an extra component, similar to that seen for S 106. In order to vi¬ 
sualize the offset caused at high-J, we have included best fit lines 
to the mid-J linear portion of the data for each obj ect in Fig. [8] 
In hig h mass protostars such as those studied by iKarska et al.l 
(l20I4l) the additional excitation is thought to arise from a com¬ 
bination of shocks and PDR emission (for low mass protostars 
such excitation is usua l ly tho ught to be dominated by shocks, 
e.g.. Ivan Kemnen et al.l 120101) . These pos sibilities will be dis¬ 
cussed in Section 16.31 In agreement with Ivan den Ancker et alJ 
(Hooi, we conclude that the majority of the CO emission ob¬ 
served towards S 106 is due to a dense PDR, particularly for the 
low to mid J lines, however it is clear that the high J CO lines 
may be indicative of a contribution by a shock. 

4.4. Summary of Diagnostics 

The simple diagnostics reveal two very similar sources, with 
subtle differences possibly driven by the different physical struc¬ 
tures of the two sources. Both show high Go, as revealed by the 
total IR fluxes, but with lower radiation fields indicated by the 
various fine-structure and molecular line ratios. Thus a two com¬ 
ponent model is required, where the PDR is comprised of two 
distinct regions: i) a hot (T ~ 400 K), dense region characterized 
by high Go (>10^) and high density (>10^ cm“^) traced by the 
various '^CO lines; and ii) a less dense region region dominating 
the emission of the atomic cooling lines characterized by mod¬ 
erate UV fields (Go ~10"^) and densities (10"^ cm“^) with a tem¬ 
perature of around 300 K. The presence of such a distinction is 
not a surprise however, as models of PDRs have inc luded density 
contras ts of this order since the early nineties (e.g.. lBurton et al.l 
Il990bl) and it has been observed in a number of objects, from 
Orion to M17 (e.g.. lStutzki & Guestenll9^lTauber et al.ll994l) . 
In addition the very high excitation CO lines for S 106 show an 
even hotter component, which may arise from a shock or a much 
hotter PDR than those considered thus far. A detailed compari¬ 
son between our sources and typical galactic PDRs will be pre¬ 
sented in Section|6T| 

5. Numerical PDR models 

The densities and UV fields found in the previous section can be 
regarded as line of sight averages as they do not take into account 
the geometry of the system or the absolute fluxes of the lines. 
In order to extract extra information from the observations, we 
used the simple analysis as the input for numerical PDR codes. 
This process should help illuminate the structure along the line 
of sight producing the emission observed. 

We use the PDR model described by IWolfire et al.l (1201 Oh . 
but modified to use either constant density or constant ther¬ 
mal p ressure. The model builds on that of iKaufman et al.l 
( 200^, which in turn comb ined the initial PDR rn odels of 
iTielens & Hollenbachl(ll985bl) and IWolfire et al.l(ll990l) . In addi¬ 
tion the model includes both the molecular hydro gen processes 
incorp orated into the Meudon PDR code (e.g., iLe Petit et al.l 
l2006l) and molecular freeze-out a nd oxygen chemistry as de¬ 
scribed bv iHollenbach et ^ (l2012l) . The model assumes a plane 
parallel slab geometry illuminated from one direction with a 


specified ID UV field and calculates the thermal balance and 
abundances of the dominant atomic and molecular species as 
well as their line emission. These models have been used to in¬ 
terpret a variety of observations of molecular and atomic species 
in a wide range of physi cal conditions including OH^, H 20 ^, 
and H 30 ^ in d iffuse gas (iHollen bach et al.ll2012£ H 2 O and HF 
in diffuse gas (ISonnentr ucker et ^ 2015h. 1C i 1 and CO in low 
UV field molecular c louds ( iBurton e t al.|l2014] : lUee et al.1l2014l) . 
H 2 in dense PDRs (ISheffer et al.l 120111) and [C ii] and [O i 
in dense PDRs and extragalactic nuclei (IKaufman et al.l UOPT 
ICroxall et al.ll20H . 


In order to get good fits to the observations of both PDRs 
constant thermal pressure models were used, with one exception 
which will be discussed shortly. In each case the final model 
shown in Fig. 0is a combination of high pressure and low pres¬ 
sure models with the parameters guided by the properties of each 
regime from the simple diagnostics. Because the simple diagnos¬ 
tics report only average properties within the beam, it is neces¬ 
sary to scale the high and low density models for each source to 
match the overall fluxes. Hence we scaled the models using the 
CO line fluxes and found good agreement with the exception of 
the highest J lines for S 106 which display additional excitation 
as discussed earlier. In order to attempt to fit these lines, a third, 
extremely high density, constant density PDR model component 
was introducec(3 The degree to which each component needs 
to be scaled can be interpreted as the Ailing factor for the high 
density model, and a rough estimate of the number of PDRs on 
the line of sight for the low density model. 

Since our models are for constant thermal pressure, both 
the temperature and density vary with depth within the cloud. 
We characterize the physical conditions where the line emission 
arises at the depth where one half of the intensity is produced. 
For IRAS 23133-(-6050, the peak CO line emission (CO 13-12) 
comes from the high pressure component atT = 217 K, n = 
1.3x10® cm“^. The [O i] 63 gm line comes from a region closer 
to the cloud surface where temperatures are higher and densities 
are lower. We And T = 2030 K, and n = 7.6x10"^ cm“^. The low 
pressure component provides about half of the observed CO 5-4 
emission at T = 28 K, and n = 9.0x10^ cm“^. It also provides 
half of the [C ii] emission at T = 530 K, n = 260 cm“^. At depths 
where these lines arise, including the high-J CO lines, grain pho¬ 
toelectric heating dominates and we And no need for additional 
and exotic heating processes. For S 106 we consider only the (a) 
model which neglects the extra, high-J component. For CO 13- 
12 we And T = 214 K, n = 1.4x10® cm“^, and [O i] 63 gm we 
And T = 2070 K, n = 7.9x10^ cm“^. For the CO 6-5 emission 
and [C II] we And T = 40 K, n = 2.7x10"^ cm^^ and T = 370 K, n 
= 1.2x10^ cm“^ respectively. Again, we And that photoelectric 
heating dominates in regions where these are formed. 

The high density component dominates the overall CO lad¬ 
der fluxes for upper levels greater than Jj, = 9, and therefore the 
Ailing factor of the high density component is well constrained 
by our flts at ~ 0.62 and ~ 0.57 for S 106 and IRAS 23133H-6050 
respectively. For the low density components factors of 4.6 and 
5.1 were found for S 106 and IRAS 23133-1-6050 respectively, 
strongly suggesting that additional low density PDRs are present 
along the line of sight. Finally, the Ailing factor of the extra 
very high density PDR component included to At the highest 
J S 106 lines was 0.006. The parameters of the various compo¬ 
nents which make up the composite models are given in Table [T] 


We refer to the two component model as S 106 (a) and the three 
component model as S 106 (b). 
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Table 7. Parameters of best-fitting numerical PDR models 


Model 

Component 

Type® 

Thermal Pressure or Density 
[K cm“^] or [cm“^] 

Go 

Filling Factor 

IRAS 23133-H6050 

1 

T 

1.5(5) 

2.0(4) 

5.1 

IRAS 23133+6050 

2 

T 

1.7(8) 

1.6(5) 

0.57 

S 106 (a) 

1 

T 

5.0(5) 

1.6(4) 

4.6 

S 106 (a) 

2 

T 

1.8(8) 

1.6(5) 

0.62 

S 106 (b) 

1 

T 

5.0(5) 

1.6(4) 

4.6 

S 106 (b) 

2 

T 

1.8(8) 

1.6(5) 

0.62 

S 106 (b) 

3 

D 

1.0(8) 

1.6(5) 

0.006 


Constant thermal pressure (T) or constant density (D). 



Fig. 9. PDR model fits to *^CO ladders, (a) Two component PDR model for S 106; (b) three component PDR model for S 106; (c) 
two component PDR fit for IRAS 23133-1-6050. In each panel the dashed line represents the low pressure component and the dot- 
dashed lines the high pressure component. In the centre panel the final, very high density PDR component for S 106 is represented 
by a triple-dot-dashed line. 


In addition to the CO line fluxes, the models also report the 
fluxes of each of the fine-structure lines in the SPIRE/PACS 
range. In Table |8] we show the ratio of the measured line fluxes 
and those predicted by the numerical models. For most lines 
the model does a good job of replicating the observed line flux 
within a factor of 2-3. As the [O i] 63 fim and [C ii] 157 /rm lines 
are commonly employed as tracers of shocks, we investigated 
whether the underprediction of the fluxes of these lines could 
be due to a contribution of a possible shock. We have estimated 
the possible contribution to these lines from shocks in S 106 by 
measuring the fraction of the total '^CO flux present in the high 
excitation component which possibly represents a shock com¬ 
pared to the dominant, cooler, components which are well fit by 
PDR models. This ratio is around 5%, showing that even in S 106 
where there is evidence for a shock in the CO lines, it is over an 
order of magnitude smaller in flux than the PDR emission in 
the CO lines. Assuming that the effect on the ionic lines would 
be of similar magnitude, shocks cannot explain the discrepancy 
between our models and observations. In fact, the differences 
in line flux estimates between different PDR models are of the 
same order of magnitude as the discrepancies we observe be¬ 


tween our numerical model fluxes and the observed line fluxes: 
therefore we conclude that the extra emission in these lines is un¬ 
likely to be due to shocks. In particular the IRAS 23133H-6050 
[C It] 157 fim line provided the most discrepant line flux ratio 
of around four, which is somewhat strange given both sources 
were fit with similar parameters and only the [C ii] in IRAS 
23I33-I-6050 shows a large deviation from the model. We sus¬ 
pect that additional low density [C ii] line emission lies along 
the line of sight to the source, or alternatively this could reflect 
complexity in the geometry of the source which could perhaps 
be resolved by higher spectral resolution measurements which 
would reveal the velocity structure. 

6. Discussion 

6.1. Comparison with Gaiactic PDRs 

While the Orion region remains the prototypical example of a 
PDR, many other examples have been observed and interpreted 
using PDR models. In Table |9] we show the properties of the 
regions studied here alongside a sample of directly comparable 
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Table 9. Summary of published properties of Galactic PDRs sorted by the spectral type of the primary exciting star. 


Object 

Spectral 

Type 

Inter-clump Phase 

Go n 

[Habing units] [cm“^] 

Dense Phase 

Go" n 

[Habing units] [cm^^] 

e 

X 10“^ 

Il58 / l63* 

Ii45 / les*” 

References 

H II Region PDRs: 

30 Dor 

03" 

7.8 X 10^ 

6.8 X 103 

- 

- 

- 

- 

- 

1, 2 

M17 

04" 

5.6 X 10^ 

3.0 X 10" 

- 

5.0 X 103 

2.3 

0.17 

0.13 

3 

M17 

04" 

- 

1.0 X 103 

- 

- 

- 

- 

- 

4 

M17 

04" 

- 

- 

1-8 X 10" 

2.0 X 103 

- 

- 

- 

5 

W49A 

05" 

3.2 X 10^ 

1.0 X 10" 

- 

- 

0.14-0.36 

0.34 

0.15 

6 

Orion Bar 

06 

4.4 X 10" 

2.0 X 10" 

- 

~103 

- 

0.11 

0.07 

7, 8, 9, 10 

Cepheus B 

OB assoc 

- 

- 

1.7 X 103 

1.0 X 10" 

- 

- 

- 

11 

S 106'' 

07-9 

1.6 X 10" 

0.6 - 3 X 10" 

1.6 X 10" 

3.1 X 10" 

1.0 

0.13 

0.15 

This paper 

S 106 

08 

- 

- 

6 103 

10" - 10" 

- 

- 

- 

12 

S 106 

07-9 

7 103 

~3 X 10" 

1-2 X 103 

1-3 X 10" 

- 

0.11 

0.17 

13 

SgrB2 

07 

5 X 103 

5 X 103 

- 

- 

- 

- 

- 

14 

IRAS 23133+6050'' 

08.5-9.5 

1.0 X 10" 

0.3-1.2 X 10" 

1.6 X 10" 

3.1 X 10" 

1.1 

0.21 

0.15 

This paper 

NGC 2024 - IRS 2b 

08-B2 

~5 X 10" 

~5 X 10" 

- 

- 

- 

- 

- 

15, 16 

S 140 

BO 

255 

2.0 X 103 

- 

2.0 X 10" 

- 

- 

- 

17 

W3 IRS 5 

B0.5 

1.0 X 10" 

1.5 X 10" 

- 

1.8 X 10" 

- 

- 

- 

18, 19 

S 125 

B0.5 

9.0 X 10^ 

1.8 X 103 

- 

- 

- 

- 

- 

20 

Reflection Nebula PDRs: 

NGC 2023 

B1.5 

1.5 X 10" 

2.0 X 10" 


- 

2.6 

0.19 

0.06 

21 

NGC 2023 

B1.5 

- 

- 

1.7 X 10" 

2.0 X 10" 

- 

- 

- 

22 

NGC 2316 

B2-3 

8.1 X 10" 

1.0 X 10" 

- 

- 

1.8 

0.42 

0.13 

23 

NGC 7023 

B3 

2.6 X 103 

4.0 X 103 

- 

- 

2.9 

0.85 

- 

24 

NGC 7023 

B3 

2.0 X 103 

7.0 X 103 

- 

- 

- 

- 

- 

25 

NGC 2245 

B7-8 

3.5 X 103 

5.0 X 103 

- 

- 

7.6 

0.75 

0.14 

23 


If different from the interclump medium radiation field. 

Observed line ratios, typically interpreted as density diagnostics. _ 

For these sources the PDR is illuminated by a cluster of stars, in the spectral type column we have indicated the hottest star detected in the literature ie.g. iHoffmeister et alJl200^ for M17 and 
iDePree etalj 19971 for W49A). 

For S 106 and IRAS 23133+6050 we quote the Interclump Go and n values indicated using the combinations of [O i] 63, 145 pm, [C ii] 158 pm and the total far-infrared flux reported in rows 3-4 
of Table|3 The dense phase (clump) values arising from the average of the high-1 CO lines and the cooling lines / FIR flux (rows 7-8 in Table O. 


References. (1) | B oschetan ( Il999h: (2) iPoglitsch et al.l (Il995h: (3) iMeixner et al.l (Il992h 

all d 1 9931) riffilHerrmannetalJ (11997b' 

C 


iTielens & Hollenbacti 


Goicoechea et al.lGZUU4I): (13 lltiiK: et al.lGZUU4I): (IblUjiannim et alj GZUUUI): (1 /(ISnaans & van JJisftoecKUiqq/l): (l»)IKj-amer et at.l(IZUU4l):( tqilMege; 

(21) ISteiman-Cameron et al.l ( ll997b : (22) ISheffer et all ( 1201 Ih : (23) lYoung Owl et al.l ( l2002b : (24) IChokshi eran ( ll988ll : (25) lJoblin et^ ( l201ol) . 


__ (4) iPerez-Beaupuits et ^ ( l201(]h : (5) ISheffer & Wolfirel (l2013h: (6) IVastel et all (l200lh: (7 ) 

(10)|Goudis|(11982 ); (11) Mookeriea et ai.l 2012); (1 2)lvan den Ancker et ai.l(l200(]h: ( 13)ISchneider et al.l(120031): 


1^85^ ; (8)mclens_rt 

2004 h (15)[^k et al.l (120031): (16 )lGiannini et alj( l200(]h: (17 )ISDaans & van D ishoeck 1 1997|); (1^ Kramer et al.l (12004 );(19) Megeath et al.l ( l200^ : (20) lAannestad & EmervI 
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Table 8. Observed / Model Line Flux Ratios. 


Line 

Wavelength 

[pm] 

IRAS 

23133-H6050 

S 106 (a) 

S 106 (b) 

[C i] 

609.14 

0.55 

0.43 

0.43 

[C i] 

370.42 

0.35 

0.29 

0.29 

[Cii] 

157.75 

4.04 

1.65 

1.64 

[O i] 

145.54 

2.43 

1.89 

1.84 

[O i] 

63.18 

0.88 

0.71 

0.71 


PDRs associated with H ii regions, as well as the properties of 
some PDRs generated by cooler star^Hl If we consider initially 
the properties of the interclump phase, our sources lie some¬ 
where between these two classes of sources, with lower overall 
radiation fields than M17, W49A or Orion - as one might ex¬ 
pect given their relatively cooler exciting stars and the fact that 
M17 and W49A are illuminated by clusters of O stars. In fact, 
the central clusters of both M17 and W49A both possess several 
examples of stars with the spectral type listed in Table |9] The 
Go values we find for our sources are actually much closer to 
those found for the reflection nebulae, particularly that of NGC 
2023. In terms of density, the interclump medium for our PDRs 
is consistent with the other examples shown in Table |3 with the 
exception of the Orion bar which presents a density around an 
order of magnitude higher. For the dense clump component the 
density deri ved for S 106 is very s imilar to that found before 
for S 106 bvISchne ider et al .l (2003h . and also close to the value 
quoted bv iMeixner et al.l(ll992h forM17. These values are typi¬ 
cal of those quoted for the clumps in a dense PDR, however the 
value of ~ 10^ cm~^ for the dense component of the Orion bar 
PDR bv lTielens et al.l (Il993h is actually very similar to the very 
high density PDR component we invoked in our numerical mod¬ 
els to attempt to fit the highest-! CO lines in S 106. This would 
suggest that the high-J component might also be present for the 
Orion bar region. Hints of this phenomenon may be visible in the 
SPIRE observations of the Orion bar presented by iHabart et al.l 
(l2010ll . where the presented RADEX fits to the *^CO lines begin 
to deviate in the mid-J regime (J ~ 12), however the region was 
not observed using the PACS spectroscopy range scan mode em¬ 
ployed for our objects, so observations with SOFIA or other FIR 
observatories would be required to confirm this hypothesis. 

Several groups have shown that there is a relationship be¬ 
tween the incident radiation held and the PDR density for a vari¬ 
ety of objects. Studies have exa mined this relationship f or a sam¬ 
ple of reflection nebula PDRs (lYoung Owl et al.ll2002l) . while a 
more recent study found a slightly different relationship for spa - 
tially resolved regions of the Carina nebula (lOberst et al.ll201 ill . 
If we consider the upper bound values for the inter-clump den¬ 
sity in Table 13 our two sources fall almost perfectly on both of 
these trends, as is shown in Figure [10] The dashed line repre¬ 
sents the expected conditions for systems in a state of equilib¬ 
rium between the pressure and temperature of the H ii region 
and PDR respectively under the assumption that the H ii re¬ 
gion is a Stromgren sphere and has an electron temperature of 
8000 K (lYoung Owl et al.ll2002t referr ed to as ‘Young-Owl et al. 


__ ling; _ 

theoretical’ in Figure [TOli. Conversely. lOberst et al.l (1201 Ih sug¬ 
gest that their best fit line represents systems which are density 


’' The optical part of such systems are usually referred to as reflection 
nebulae, as the primary optical emission from such nebulae consists of 
reflected starlight rather than emission lines as one would expect from 
H II regions. 


bounded as opposed to the ionization bounded RNe, and hence 
they found a different relationship. We have also included the 
H II region PDR sources from Table [9] while the R Ne sources 
in Table 13 form part of the lYoung Owl et ^ (l2002h sample and 
are represented by the black diamonds in Figure[T3 Our sources 
are consistent with the H ii region PDRs with late O type ex¬ 
citing stars and the RNe formed by early B stars. Amongst the 
H II region PDRs with hotter exciting stars there seems to be a 
greater spread in properties away from the trends, with W49A 
and W3 IRS 5 in particular having either higher Go or lower hh 
by at least an order of magnitude than predicted by either rela¬ 
tionship. We have also included the properties of the dense phase 
for several sources, denoted as ‘dc’ in Figure[T0| For our sources 
and NGC 2023, the properties of the dense phases fall very near 
either of the relationships discussed, the dense component of 
M17 ho wever falls well abo ve even the matter bounded trend 
given by lOberst et al.l (1201 111 by several orders of magnitude, 
this may be due to the dense clumps being independently gravi¬ 
tationally bound or transient objects - that would therefore have 
much higher densities than would be expected for pure pressure 
equilibrium between the clumps and their surroundings. Overall 
it seems that the simpler sources with a dominant exciting star 
(e.g. IRAS 23I33-H6050, S 106, NGC 2023, S 125, S 140) agree 
well with the trends, while the much more complex sources (e.g. 
W49A, M17, W3) have much more variance. However one of 
the most complex sources, 30 Dor, is consistent with both pre¬ 
dictions, so the possible trend between simple/complex sources 
could easily be a manifestation of small number statistics. 

It is expected that the flux ratio of the 158 /rm [C ii] and the 
63 jum [O i] lines should decrease with increasing Go and n as the 
[O i] line becomes the dominant source of cooling. We see this 
progression between our two sources, with IRAS 23133-(-6050 
having a cooler ionizing star and lower Go and a correspond¬ 
ingly higher / j-Q jj ratio than S 106 by around a factor of 
two. The two values we And for this ratio fit nicely into the pro¬ 
gression seen in Table |3 with lower values for the H ii regions, 
(with the puzzling exception of W49A although this may be due 
to the large physical size of W49A in comparison to the other 
objects) and higher values for the RNe PDRs. As with the com¬ 
parisons between the radiation held and densities of the regions, 
the closest Galactic analog of IRAS 21333-1-6050 we measured 
is actually the reflection nebula NGC 2023. This likely arises 
because the 08.5-9.5 star illuminating IRAS 21333-1-6050 has a 
surface temperature relatively close to that of the B1.5 central 
star of NGC 2023. On the other hand, S 106 is closer to the H ii 
region PDRs listed, and actually shows the same line ratio as the 
Orion bar region. 

In context of other Galactic PDRs then, on the basis of the 
spectral types of the exciting stars alone the PDRs studied here 
might be expected to have properties between those of the mas¬ 
sive star forming regions and clusters and the RNe PDRs, and 
this is indeed what we And. With both regions having properties 
similar to the PDRs generated by the hottest B stars, especially 
IRAS 23133-(-6050, which is very similar to the PDR of the RNe 
NGC 2023. The general consistency between our two objects 
can also be explained using similar logic. Both sources have sim¬ 
ilar H II region electron densities and sizes, and similar types 
of exciting stars, around 08-9, although from the literature pre¬ 
sented in section|3it seems likely that the star illuminating S 106 
is slightly hotter by about one subtype. Despite this, both sources 
have very different morphologies, with S 106 having a very com¬ 
plex dual lobed H ii region and a small stellar disk region sur¬ 
rounded by a large extended PDR, while IRAS 23133-1-6050 is 
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log G(j 


Fig. 10, Correlat i ons be twee n Gn and de n sity w p found by 
lYcmng Owl et al.l (l2002h and ObersteLaf (1201 Ih . Reflection 


nebula points from lYoung Owl et al.l (120021), Galactic PDRs as 


stated in Tablel9] IRAS 23133+6050 and S 106 fit into the trend 
between the larger H ii regions with hotter exciting stars and the 
RNe PDRs with B type exciting stars. Dashed and solid black 
lines: theoretical prediction for the relationship between radia¬ 
tion field and density for a pressure b ound H ii region and the 
best fit to the lYoung Owl et al.l (l2002h da ta. Green dot - dashe d 
line: best fit to the Carina nebula data of lOberst et al.l (1201 ll) . 
For a selection of sources we have also included the points rep¬ 
resenting the dense PDR component (denoted ‘dc’) as opposed 
to the default values representing the interclump medium (de¬ 
noted ‘ic’). 


still an ultra-compact H ii region yet to break out of its natal 
molecular cloud. This dichotomy poses an interesting problem 
and might constrain models for the evolution and dynamics of 
H n regions and PDRs. 


6.2. S 106 

S 106 in particular has a rich vein of literature and prior stud- 
ies, of which a few discuss its P DR speciflcally (for ex ample 
Ivan den Ancker et akl (l2000h and ISchneider et alJ (l2003l) listed 
in Table |3- We discussed these studies in terms of comparing 
the measured surface brightnesses for the various PDR cooling 
lines in Section [3.3.1l here we discuss the overall results of those 
studies and how they compare to our own. In general there is 
broad agreement between the literature findings and our results 
for S 106, in that there is a clumpy PDR consisting of dense 
clumps irradiated by a strong UV field producing the high-J CO 
lines, along with a less dense PDR phase permeating the re¬ 
gion. We have extended these findings by presenting the full CO 
ladder available to the PACS and SPIRE Herschel instruments, 
which have revealed the presence of a further, even higher den¬ 
sity, component. The key difference between the prior studies 
and our results here is that we are probing only the innermost 


region of S 106 (a circular aperture of radius 21.5"), whereas 
Ivan den Ancker et al.l (l2000l) used ISO-LWS spectra of the cen¬ 
tral region with a beam som e three to five times la r ger tha n our 
Herschel observations. The Ivan den Ancker et al.l (l2000l) con¬ 
cluded that the S 106 PDR was of very high density and ex¬ 
perienced an extremely intense UV radiation held (Gq^ 10^) 
despite measuring a similar set of the common PDR cooling 
lines as we presented earlier from the Herschel observations and 
found evidence of a more moderate environmen t. The key differ¬ 
ence arises because lvan den Ancker et ^(l2000l) also considered 
ISO-SWS observations including the [Si ii] 34.8 /rm line, the 
strength of which could only be explained by invoking stronger 
radiation fields and higher densities. _ 

A subsequent study by ISchneider et al.l (l2003h spatially re¬ 
solved the larger scale structure of S 106 PDR cooling line emis¬ 
sion with pixel scales of 20-40"using the FIF I instrument on 
the KA O, and also including ISO-FWS spectra. ISchneider et al.l 
(l2003h found density variations within the PDR in a similar man¬ 
ner to those described in this work, and very similar parame¬ 
ters for the two phases were recovered. The only slight differ¬ 
ence arises in the density of the dense compon ent, in which we 
found a value slightly above the range quoted bv ISchneider et al.l 
(l2003h . We also And a value for the average UV radiation held 
experienced by t he inte rclump medium of Go = 1.6 x 10^ while 
ISchneider et al.l (l2003h give a lower limit of Go ^ 10^. The fact 
that we recover very similar values for the density and incident 
radiation held for the interclump medium in the region imme¬ 
diately adjacent to the star as was found by previous studies as 
averages for the entire region can then be easily explained by in¬ 
spection of the maps presented bv ISchneider et al] (l2003l) which 
show that the emission from S 106 is very strongly peaked near 
the exciting star in the center of the region. By using the full 
Herschel CO ladder we have extended this picture to include a 
further excitation component, which may represent a further, yet 
higher density, PDR component or a shock - a conundrum we 
will discuss in the following subsection. 

6.3. The diagnostic vaiue of the CO iadder 

We have used the observed CO fluxes in order to constrain the 
PDR properties of our sources. The shape of the CO ladder (as 
well as the curves in the rotation diagrams) provides the number 
of components with distinct physical conditions present in each 
source. In combination with the simple diagnostics which use the 
standard PDR lines, the ladder was then used to find the temper¬ 
ature of each component using RADEX models. Subsequently we 
used the CO fluxes to guide the application of numerical PDR 
models - leading to a reasonable fit to all of the other observ¬ 
ables. However, there are several complexities to this picture, 
firstly: to what extent do the results of the rotation diagram anal¬ 
ysis make sense in terms of our more complex models? It is 
not clear that the level populations are in FTE while the lower 
J lines may be optically thick and the temperature distribution 
is not constant, so it might be expected that the results we And 
for the rotation diagrams compare poorly to the more detailed 
approaches. Secondly: we have detected the presence of an extra 
excitation component in the high-J S 106 CO lines which could 
be ascribed to either very-high density/thermal pressure PDR, or 
a strong shock. 

The rotation diagram fits to both sources revealed very sim¬ 
ilar rotation temperatures of around 190 K (188 and 190 K for 
IRAS 23133+6050 and S 106 respectively) for the mid-J '^CO 
lines. In Section l4.3.1l we compared these temperatures to those 
found using RADEX fits to the respective ^^CO ladders. Using a 
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density of 10^ ® as indicated by the diagnostic diagrams 
gave temperatures of around 100 K warmer. We also found rea¬ 
sonable fits to the CO ladders using RADEX with higher densities 
and a temperature of 200 K. In light of the results of our nu¬ 
merical models, these results are actually reasonable estimates. 
For example, for IRAS 23I33H-6050 we found that the peak CO 
line emission emanated from a region ofT = 217Kata den¬ 
sity roughly three times higher than that predicted by the simple 
diagnostic diagrams. This concurrency suggests that the skepti¬ 
cism with which we referred to the rotation diagram figures was 
perhaps unwarranted; in fact, for these data one could use the ro¬ 
tational temperature found using the rotation diagrams as a guide 
for using RADEX models to find the PDR density in the vicinity 
of the CO molecules. 

In a more general sense the CO ladders we find are typi¬ 
cal, and are routinely found in extragalactic contexts. This is be¬ 
cause any observation of a star-forming galaxy will inevitably in¬ 
clude star forming complexes like S 106 and IRAS 23133-1-6050. 
However, such observations may well include regions dom¬ 
inated by other physical processes, including analogues of 
PDRs dominated instead by X-rays or cosmic rays (XDRs or 
CDRs), PDRs with enhanced mechanical h eating (mPDRs) or 
an environment dominated by shocks (e.g., van der Werf et alJ 
l2010t iRangwala et al.l 1201 it iRosenberg et ah 2014 and refer¬ 
ences therein). Each of these processes has some effect on both 
the relative strengths of the CO ladder and the fine-structure 
lines. The challenge in interpreting such observations then is to 
tease apart the different effects and attempt to create a coherent 
model of the environment. In contrast our sources are relatively 
simple and do not seem to be dominated by the wide variety of 
other physical processes that can be important for star-forming 
galaxies or AGN. Therefore, the extreme conditions we infer in 
our sample represent a baseline for the more complex sources, in 
that such conditions would be expected in the PDRs associated 
with massive star formation in these extragalactic sources. 

There appear to be two plausible explanations for the extra 
heating in the very high-J (J > 20) CO lines seen in S 106; shocks 
or a very high density PDR. Previo us work on S 106 with r espec t 
to shocks is somewhat conflicting. Ivan den Ancker et al.l (l2000li 
found that the IR emission from S 106 was dominated by PDR 
emission on the basis of the [S i] 25.2 //m / [Si ii] 34.8 /tm line 
ratio, which is usually thought to be rep resentative of shocks. 
However, the Ivan den Ancker et al.l (l2000l) study used ISO-SWS 
observations with a comparatively large beam - possibly imply¬ 
ing that the environment is dominated by PDR material, rather 
th an exclusively PDR material. A later spatially resolved study 
by iNoel et akl (l2005l) found evidence of thermalization in near- 
IR H 2 line ratios, which could be evidence of shocks in the S 106 
region, which would be included in the beam of our observa¬ 
tions. Such rat ios could also be pr oduced by a very high den- 
sity PDR (e.g.. lBurton et al.lll990ah . although iNoel et alT (120051) 
conclude that this cannot be the case for one of the thermal- 
ized regions, as it lacks detections of Bry, PAHs and contin¬ 
uum emission. As we noted earlier, sho cks are common in star 
formation regions of most masses (e.g.. lSantangelo et al.ll2013t 
iKarska et alfcoidh . The bipolar morphology of S 106 lends fur¬ 
ther credence to the idea that shocks are present in S 106 and are 
contributing to the extra excitation obseryed in the CO ladder. 

In contrast, the eyidence supporting the alternatiye yiew - 
that the extra excitation arises from the high density PDR itself 
- is somewhat limited and circumstantial. If taken at face yalue, 
the center panel of Fig.|3and the negligible differences in terms 
of the fine-structure line fluxes between the two S 106 models 
in might appear to support the idea that the emission arises from 


a yery high density PDR. Howeyer, we found that a good fit to 
the high J lines could be only reached using an additional PDR 
model with a density of 10^ cm“^ and Go of 1.6 x 10^. Such den¬ 
sities could perhaps be plausible in two circumstances. Firstly, 
the small circumstellar disk around IRS4 in the core of S 106; 
secondly, small dense clumps interspersed within the main PDR. 
Howeyer, both of these scenarios could be problematic. The first 
possibility would require a much higher yalues of Go than used 
in our model by yirtue of being much closer to the primary star, 
howeyer the second circumstance would require small clumps of 
PDR material with densities seyeral orders of magnitude beyond 
any preyiously obseryed. Clearly further studies will be required 
to decisiyely settle this issue. Neyertheless, the g eneral shape 
of the C O ladder fluxes is yery similar to Fig. 5 of iBurton et al.l 
(1199^, in which at high densities and radiation fields, the high- 
J CO lines can be excited in a thin layer on the surface of PDRs 
at temperatures of seyeral hundred K. The transition between the 
two regimes is not at precisely the same howeyer the shape 
of the curye is almost identical. 


7. Conclusions and Summary 

We haye obtained full Herschel -PACS and -SPIRE spectra of 
the star forming regions S 106 and IRAS 23I33H-6050 and an¬ 
alyzed them in terms of the emission from their PDRs. Simple 
PDR diagnostic diagrams were employed to And the physical 
conditions and hints of the physical structure of the two PDRs. 
These results were tested using a detailed numerical model of 
PDR physics. The study reyeals that the PDRs around S 106 and 
IRAS 23133-1-6050 are remarkably similar in some respects, de¬ 
spite being dramatically different types of objects - S 106 is a 
large, complex star formation region, with many components, 
while IRAS 23133-1-6050 is a relatiyely simple cometary ultra¬ 
compact H II region. Despite these differences, both objects were 
found to haye a similar range of physical conditions and struc¬ 
tures, in common with many other PDR obseryations and mod¬ 
els. The structure of each region, implied by our obseryations, 
is that of a classical PDR with intermediate properties (Go=10"^, 
«=10"^ cm“^) and with embedded clumps of much denser ma¬ 
terial (Go=10^, n-lQ^ cm“^). We employed numerical models 
which were fitted to the *^CO ladder fluxes, initially to deter¬ 
mine the temperature of the PDR enyironments from the peak 
of the CO ladder, as well as to inyestigate whether the condi¬ 
tions indicated by the simple diagnostics were consistent with 
those indicated by a full numerical treatment. During this pro¬ 
cess it was noticed that the S 106 CO ladder possesses an extra 
excitation component at the high J lines. This component is only 
reconcilable with PDR models upon inyoking extreme densities, 
so we tentatiyely ascribe shocks as a probable cause. Ultimately 
though, both PDRs studied possess yery similar properties to 
those of the Orion complex, which is surprising as Orion and 
the two objects studied here are in yery different eyolutionary 
stages. The youngest object, IRAS 23133H-6050, is still embed¬ 
ded in its natal cloud, while in Orion the cluster of exciting stars 
haye emerged from this phase and are steadily eroding their en- 
yironment. The consistency of properties on yarious PDRs be¬ 
tween objects and in time underlines their importance as a key 
phase of the ISM. 
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